New type of parametrizations for parton distributions 

A. Yu. lUarionov 
Dipartimento di Fisica, Universita' di Trento, Ital'^ 

A. V. Kotikov, S. S. Parzycki, and D. V. Peshekhonov 

Joint Institute for Nuclear Research, Russia 

(Dated: May 21, 2010) 

New type of parametrization for parton distribution functions, based on an exact their Q - 
evolution at large and small x values, is constructed for valence quarks. It preserves exactly the low 
X and large x asymptotics of the solution of the DGLAP equation and obeys the Gross-Llewellyn- 
Srnith sum rule. 
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I. 



INTRODUCTION 



The parton distribution functions (PDFs) which con- 
tribute to the LHC processes, and the PDFs fitted at 
HERA and fixed target experiments are defined at some- 
what different (cc, Q^) ranges Q (see, for example. Fig. U] 
taken from 0]). Therefore, a direct application of the 
modern PDF sets |l!] to the LHC processes is not well 
justified. 

This problem is really important, because the larger 
uncertainties for many processes at the LHC originate 
mainly from our restricted knowledge of the parton dis- 
tributions (see, for example, the recent paper 3] and ref- 
erences therein). 

In the present paper we propose an idea to over- 
come this problem. Our solution consists of the two 
basic steps. Firstly, we find asymptotics of solu- 
tions of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) equations [^ for the parton densities at low 
and large x values of the Bjorken variable and, at the 
next step, we approach a combination of the two solu- 
tions for the full range of x. 

In a sense, this is not a new idea. A similar approach 
had been given by the Spanish group (see book @ and 
references therein) about 40 years ago. However, in the 
present paper the parametrization will be constructed 
in a rather different way. In particular, it includes an 
important subasymptotic term which is fixed exactly by 
the sum rules. 

The results of the present paper are restricted by the 
leading order (LO) of the perturbation expansion, what is 
reasonable Q, since for many processes at the LHC the 
next-to-leading-order (NLO) corrections are not known 
so far. Moreover, in the present study we limit ourselves 
to the valence quarks, whose evolution does not contain 
contributions from the gluons. 
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LHC parton kinematics 



■ x,2 = (M/14TeV)exp(±y) 
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FIG. 1. The (x,Q'^) ranges of the PDFs contributed to the 
LHC processes and ones fitted at HERA and fixed target ex- 
periments (from 2]). 



II. 



A SHORT THEORETICAL INPUT 



* lillario@science.umtn.itl 

^ Several recent PDF fits, as well as the references to the previous 
studies, can be found in [y. 



In this section we briefly present the theoretical part 
of our analysis. The reader is referred to [jI for more 
details. 

The deep- inelastic scattering (DIS) ^ + iV -^ T + AT, 
where / and TSS are the incoming lepton and nucleon, and 



I' is the outgoing lepton, in the one of basic processes 
to study of the nucleon structure. The DIS cross-section 
can be spht to the lepton L'^^ and hadron F'^'^ parts 



da - L^"'F^"' 



(1) 



The lepton part L'^^ is evaluated exactly, while the 
hadron one, F'^'^ , can be presented in the following form 



piiu 



qf'q" 



Fi{x,Q^) 



^^+^^A(p^+^-4i')f^M'' 



«e,,i.a/3PV^-F3(a;,0^) + 



(2) 



where the symbol . . . stands for the parts, which depend 
on the nucleon spin. The functions Fk{x,Q^) (hereafter 
k — 1,2,3) are the DIS structure functions (SFs) and q 
and p are the photon and parton momenta. Moreover, 
the two variables 



-\2 2 ^ n ^ 

) = — 9 > , X = 



2ipq) 



(3) 



determine the basic properties of the DIS process. Here, 
Q^ is the "mass" of the virtual photon and/or Z/W bo- 
son, and the Bjorken variable a; (0 < a; < 1) is the part of 
the hadron momentum carried by the scattering parton 
(quark or gluon). 



Mellin transform 



(a — qi), antiquarks (a — qi) with {i = 1, ■ ■ • ,6), and 
gluons {a — G), i.e. 



Aa{n,^i') = fain,^i')= dxx'^-^fa{x,^i'). (6) 

The coefficient functions C^{n,Q'^/^'^) are represented 
by 



C-,in,Q'/^,')^ dxx^-'Cl{x,Q'|^i') (7) 
Jo 

and responsible for the relationship between SFs and 
PDFs. Indeed, in the a;-space the relation ([S]) is replaced 
by 

Fk{x,Q^)= J2 Cl{x,Q^/^i'')®fa{x,^i^), (8) 

a—q,q,G 

where (^ denotes the MeUin convohition 



Jx y \y 



(9) 



Using Eqs. ([5]) and ([8]), one can fit the shapes of 
PDFs /a(a;, ^^), which are process-independent and use 
them later on for another processes. Indeed, the cross- 
sections of the hadron-hadron processes are proportional 
to the parton luminosities fafi{x,p?) (see, for example, 
the recent paper 8] and references therein), which are 
the Mellin convolutions of the two PDFs fa{x,fM^) and 
fb{x,fi^): 



fa.b{x,tx^) = Ia[x,^l^) (g) fb{x,^i'^). 



(10) 



The Mellin transform diagonalizes the Q^ evolution of 
the parton densities. In other words, the Q^ evolution of 
the Mellin moment with certain value n does not depend 
on the moment with another value n' . 

The Mellin moments Mk{n, Q^) of the SF Fk{x, Q^) 



Mk{n,Q^)= / dxx''~'Fk{x,Q^) 



(4) 



can be represented as the sum 



Mfe(n,Q2) 



a—q,q,G 



Cfe"(n,QVM')^a(n,/i'), (5) 



where C^{n,Q^/ii^) are the coefficient functions and 
Aa{n,fi'^) =< N\0'^_^ nJ^ > ^^^ the matrix elements 
of the Wilson operators O" , which in turn are pro- 

cess independent. 

Phenomenologically, the matrix elements Aa{n, /i^) are 
equal to the Mellin moments of the PDFs fa{x,fi^), 
where /a(a;,/i^) are the parton distributions of quarks 



^ All parton densities are multiplied by x, i.e. in the LO the struc- 
ture functions are some combinations of the parton densities. 



B. Quark distribution functions 

The distributions of the u and d quarks contain the 
valence and the sea parts 

fqi ^ fu = fu '^ fu J fq2— fd = Jd^Id- (H) 

The distributions of the other quark flavors and of all the 
antiquarks contain the sea parts only: 

f,.=fqp (J = 3,4,5,6), /g. =/| (z = l,---,6). 

(12) 

It is useful to define the combinations of quark densi- 
ties, the valence part /y, the sea one fs and the singlet 
one fsi i 

6 
fv ^ fu + fd ^ fs ^22 (-^1 + -^1) ' 



/s/-^(/,. + /,-.) = /v + /s. 



(13) 



j=l 



^ Here we consider all quark flavors. Really, heavy quarks factorize 
out when \/Q^ becomes less then their masses, and we should 
exclude them from the Q^-region. 



Because the PDFs, which contribute to the structure 
functions, are accompanied by some numerical factors, 
there are also nonsinglet parts 



fA.,^{U+U)~{fg,+h) 



(14) 



which contain difference of densities of quarks and anti- 
quarks with different values of charges. 

As an example, we consider the electron-proton scat- 
tering, where the corresponding SF has the form 

6 

Ff (x,g2) = ^e,2 {U{x,Q^) + f,^{x,Q')) . (15) 

4=1 

In the four-quark case (when 6 and t quarks are separated 
out), as in [9j, we will have 



F2'i^, Q') = ^ fsiix, Q') + i/A(x, g2) ., 
where 

fA= J2 iUi^^Q') + Ui^^Q')) 

q-i—u.c 
qi—d,s 

C. DGLAP equation 

The PDFs obey the DGLAP equation [i| 
d 



(16) 



(17) 



dln/i' 



fa{x,H^) ^J2^''''{x)(E)fbix,^I^), (18) 



where a,b — NS, SI, G and ^ab{x) are the so-called split- 
ting functions. 

Anomalous dimensions ^ab{n) of the twist-two Wilson 
operators Of, ,, in the brackets h are the Mellin trans- 
forms of the corresponding splitting functions 



lab{n) = dx x" %b{x) . 
Jo 



(19) 



In the Mellin moment space, the DGLAP equation be- 
comes to be the standard renormalization-group equation 



fa{n,fi^) = E ^ab{n)fb{n,fi^). (20) 



dln/i^ 



b=NS,SI,G 



Below we will study the properties of the valence part 
only. Consideration of the other quark densities will be 
discussed separately. 



III. LOW AND LARGE x ASYMPTOTICS 

The large x asymptotics of the valence quark density 
has the following form [l^, llll 

fv{x,Q') ^ By{s){l-x)^-^-^\ (21) 



where 

f\n(Q'^/A'^)\ - - 16 

' ^ ^" (,ln(QVA2)j ' ^^(^) = ^^(0) + ^^^ ' '^^ = 3^ ' 



Bvis) = Bv{0) 



rii + i3vis)) 



Pv 



dv [ IE - 2 ) s . 



(22) 



Bv{0) and /?v(0) are free parameters. Here /3o = 11 — 
2//3 is the first term of the QCD /3-function, / is the 
number of active quarks and je is the Euler constant. 
The constant PyjO) can be estimated from the quark 
counting rules [l2| as 



/3i/(0)~3. 



(23) 



Eqs. ([21]) and (|22|) demonstrate the fall of the parton 
densities at large x values when Q^ increases. 

At small- a: values the valence part has the following 
asymptotics [H 0, [11 



where 



fv{x) -^ Avis)x^ 



Ay(s) = Av.(0)e-'^"^(i~^^)% 



(24) 



(25) 



dNsin) 



16 



*(n + l) +7B 



2n{n + 1) 



Xy and Ay{0) are free parameters and $(n+ 1) is Euler 
^-function. 

From the Regge calculus the constant Ay ~ 0.3 -f- 0.5. 
Moreover, the Q^ evolution of this parton density shows 
that \v should be Q^ independent [15]. 

In Eq. (|25p the "anomalous dimension" dNsin) is rep- 
resented in the form which is useful for the non-integer n 
values. Usually, the elements of the coefficient functions 
and anomalous dimensions are expressed as the combi- 
nations of the nested sums [la], which can be expanded, 
however, to the non- integer n values according to |17|. 



IV. PARAMETRIZATION 

The valence quark part fv{x,Q'^) can be represented 
in the following form|j 



fv{x,Q^)=x^^{l-x)^^'''^ 



Avis){l-x) +Bvis)x 



(26) 



* Following to [isll . it is possible to add to the parametrization 
1261 1 an additional polynomial term (l + 5Zfe=i '^V,k^^) with cer- 
tain constants ay^k- This term fixes properly the PDF shape 
and may improve essentially an agreement with the correspond- 
ing experimental data. In the valence sector, considered here, 
we found a good agreement between II26I I and the data without 
such term. However, in the singlet part, where there is a strong 
correlation between the see quark and gluon densities, similar 
terms may be useful. 



which is constructed as a combination of the small x 
and large x asymptotics, the last term is equal to Ay{s) 
at small x and to Bv{s) at large x values. The Q^- 
dependence of the parameters in (j26p is given by Eqs. 
and (1251). 



A. Gross-Llewellyn-Smith sum rule 

The additional relation between the parameters in (j^ 
stems from the LO Gross-Llewellyn-Smith sum rule ^1<^ 







where the last term in the brackets is subasymptotic in 
both the small and the large x values. These types of the 
subasymptotic reduction, ~ ^/x at small x and ^ (1 — x) 
at large x values, have been discussed in [5]. 

Now, the sum rule ^7} can be satisfied at any Q^ 
values and determines completely the new term. 



Qi 



T{\v)T{l + l3v{s)) 
T{Xv+2 + l3vis)) 
r(Av. + l/2)r(2 + /3y(g)) 
r(Ay + 5/2 + /?y(s)) 



XvBv{s) + {1 + pv{s))Av{s 
Dy{s). (31) 



^ dx 



fv{x,Q^)^Qv , gy = 3. 



(27) 



which simply informs about the number of the valence 
quarks in nucleon. 

As long as the Eq. (1261) is just a parametrization, an 
attempt to apply the sum rule for, e.g., k different values 
of Q^ produces k additional relations that is, of course, 
unacceptable. The sum rule can be applied only at one 
point s = sc for some critical value of Q^. For the other 
values of Q^ , the sum rule holds only approximately and 
one can estimate the deviation from the exact sum rule. 

So, we have the following relation 



Qv 



T{Xv)r{l+l3v{sc)) 



AyBy(0)e-P^^ 



r(Ay + 2 + /3y(sc)) 

+ (1 + /3v(sc))Ai.(0)e-'^"^(i-^^)^^ 



(28) 



It is possible to choose any "middle" value of s. How- 
ever, if sc — 0, i.e. Qq — Qg, then Eq. (^51) is consider- 
ably simplified 



Qv 



r(Ay)r(i + /3v(o)) 

r(Ay + 2 + /?y(0)) 



XvBviO)- 



■(l+/3y(0))Ay(0) 

(29) 



V. RESULTS 

To obtain the parameters of our parameterization (|30p 
for the valence part, we compare it numerically with the 
results of several Gluck-Reya-Vogt (GRV) setsjp, [12] 
and Gluck-Jimenez-Delgado-Reya (GJR) one [2J|. The 
choice of the GRV/GJR evolution is related to future in- 
vestigations of the gluon and sea quark densities; the 
small X asymptotics obtained in [24\ are conceptually 
close to the GRV/GJR approach. 

With the parameters, reported in the Table |T] our pa- 
rameterization (pO| and the GRV/GJR sets are in good 
agreement for all x values in the very broad Q^ range, 
0.25 GeV2 < Q2 < 194 QeV^. 

Then, the parameterization pop describes the experi- 
mental data as well as the GRV/GJR sets and has the 
analytic form containing the exact asymptotics (|21l) and 
([M)l of the DGLAP evolutions at small and large x values. 
Thus, it can be applied with good warranty in any other 
(x, Q^) region, where the experimental data are still not 
available. So, it should be applicable for the LHC range 
of X and Q^ values. 

Note, that all GRV/GJR sets themselves give quite 
close results for the valence quarks, i.e. these results 
should not change significantly when the new data will 
appear, for example, from the LHC experiments. 



B. Subasymptotic term 



VI. CONCLUSIONS 



In the present paper we choose another possibility to 
take exactly the sum rule (l?7|) into account. We intro- 
duce an additional term ~ Dv{s) to our parameteriza- 
tion (j26p . which can be written as 

fv{x,Q^)^x^^{l~xy^^'^ (30) 

X Avis)il-x)+Bvis)x + Dvis)^{l-x) , 



In this work, we investigated the low x and large 
X asymptotics of the valence quark density, and ob- 
tained the corresponding parametrization pop . This 
parametrization obeys explicitly the Gross-Llewellyn- 
Smith sum rule. It has been compared with various 
GRV/GJR sets ^^^ in order to fit the initial values of 
all the Q^-dependent parameters. It has been performed 
accurately, because numerically the form (150]) and the 
GRV/GJR sets have very similar shapes for all x and Q^ 
values. 



^ Above LO, the Gross-Llewellyn-Smith sum rule |T9'| is defined 
as the integral of the structure function _F3 and contains the 
perturbative (~ as) and the power corrections in its r.h.s. (see, 
for example, |20l| and references therein). 



We used Q^ evolution at fixed / value. The contributions of the 
heavy-quark thresholds are negligible | 25| . 





/ 


Q.(Afz) 


Ql (GeV^) 


Av{Q) 


Bv{Q) 


Av 


/3v(0) 


GRV(1991) [21J 
GRV(1998) [22] 
GJR(2008) [23] 


4 
3 
3 


0.1282 
0.1250 
0.1263 


0.25 
0.26 
0.30 


0.61 
0.93 
1.07 


80.62 
64.80 
85.60 


0.317 
0.333 
0.390 


2.870 
2.812 
2.945 



TABLE I. The parameters of the valent quark distribution of Eq. (|30l) fitted on several GRV/GJR sets. 
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FIG. 2. The valence quark density at Qq- 



At the next step, we plan to add to our analysis the 
nonsinglet and sea quark densities as well as the gluon 
distribution, and apply the obtained results to the anal- 
ysis of several LHC processes. Moreover, a comparison 
of these results with the predictions of other sets [l| of 
the parton densities will be done. We also plan to con- 
sider the PDFs in nuclei. Therefore, the EMC effect p6|, 
which is very important in the high-energy regime, will 
be added to the consideration. 
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